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Abstract 

A novel method is developed for extending the Green-Naghdi (GN) shallow-water 
model equation to the general system which incorporates the arbitrary higher-order dis¬ 
persive effects. As an illustrative example, we derive a model equation which is accurate 
to the fourth power of the shallowness parameter while preserving the full nonlinearity of 
the GN equation, and obtain its solitary wave solutions by means of a singular pertur¬ 
bation analysis. We show that the extended GN equations have the same Hamiltonian 
structure as that of the GN equation. We also demonstrate that Zakharov’s Hamiltonian 
formulation of surface gravity waves is equivalent to that of the extended GN system by 
rewriting the former system in terms of the momentum density instead of the velocity 
potential at the free surface. 
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1. Introduction 


The Green-Naghdi (GN) equation which is also known as the Serre or Su-Gardner 
equations models the fully nonlinear and weakly dispersive surface gravity waves on fluid 
of hnite depth. See Serre [1], Su and Gardner [2] and Green and Naghdi [3]. Although the 
GN equation approximates the Euler equations for the irrotational flows, it exhibits several 
remarkable features. In particular, it has a Hamiltonian formulation which provides a 
unified framework in exploiting the mathematical structure of various model equations 
such as the Boussinesq, Korteweg-de Vries (KdV) and Gamassa-Holm (GH) equations 
(Gamassa and Holm [4] and Gamassa et al [5]). A large number of works have been devoted 
to the studies of the GN equation from both analytical and numerical points of view. A 
review artile by Barthelemy [6] describes the derivation of the GN equation as well as a 
method for improving the dispersive effect. Furthermore, the improved model equations 
are tested agianst experiments. The recent article by Bonneton et al [7] reviews the high- 
order numerical methods for the GN equation and the numerical results in comparison 
with breaking random wave propagation experiments. The following two monographs are 
concerned with the derivation and mathematical properties of the GN and other water 
wave equations: Gonstantin [8] provides an overwiew of some main results and recent 
developments in nonlinear water waves including breaking waves and tsunamis. Lannes 
[9] addresses the derivation of various asymptotic model equations and their mathematical 
analysis which is mainly devoted to the well-posedness of the model equations. 

The GN equation incorporates the dispersion of order <5^, where 6 = ho/1 is the shal¬ 
lowness parameter {Hq: mean depth of the fluid, 1: typical length scale of the wave). To 
improve the dispersion characteristics, various attempts have been made to extend the 
GN equation. Among them, the model equations have been derived which include the 
dispersive terms of order <5^ (Kirby [10], Madson and Schaffer [11, 12] and Gobbi et al 
[13]). Numerical computations have been performed for these equations to examine the 
wave prohles and the amplitude-velocity relations as well as the effect of dispersion on the 
wave characteristics. Note, however that whether the proposed higher-order dispersive 
model equations permit the Hamiltonian formulations has not been discussed so far. In 
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this paper, we extend the GN equation which is accurate to the dispersive terms of order 
while preserving the full nonlinearlity, where n is an arbitrary positive integer. The 
case n = 1 corresponds to the GN equation. We show that the extended model equations 
have the same Hamiltonian structure as that of the GN equation. 

We consider the two-dimensional irrotational flow of an incompressible and inviscid 
fluid of uniform depth. The effect of surface tension is neglected for the sake of simplicity. 
The governing equation of the water wave problem is given in terms of the dimensionless 
variables by 


-|- (pyy = 0, —oo < X < oo, —l<y<eri, 

( 1 , 1 ) 

m + = ^(/y, y = 

( 1 , 2 ) 

+ + = y = e7], 

( 1 , 3 ) 

(py = 0, y= -1- 

( 1 , 4 ) 


Here, 0 = (f){x,y,t) is the velocity potential, y = ri{x,t) is the prohle of the free surface, 
and the subscripts x, y, and t appended to 0 and rj denote partial differentiations. The 
dimensional quantities, with tildes, are related to the corresponding dimensionless ones 
by the relations x = lx,y = h^y^t = {l/co)t,fj = arj and (p = {gla/co)(f), where a and Cq 
are characteristic scales of the amplitude and velocity of the wave, respectively, and g is 
the acceleration due to the gravity. Note that Cq = y/gho is the long wave phase velocity. 
In the problem under consideration, one can choose the two independent dimensionless 
parameters, e = a/ho and S = ho/1- The former parameter characterizes the magnitude 
of nonlinearity whereas the latter characterizes the dispersion or shallowness. 

In Section 2, we provide a recipe for deriving the model equations. See, for instance 
Matsuno [14] as for an analogous method which develops a procedure for obtaining the full 
dispersion model equations of the water wave problem. After completing the construction 
of the extended GN system, we derive as an example a model equation which is accurate 
to order 5^. In Section 3, we show that the extended GN equations can be formulated 
in a Hamiltonian form by employing an appropriate Lie-Poisson bracket. At the same 
time, we demonstrate that the extended equations are equivalent to Zakharov’s equations 
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of motion for surface gravity waves. In Section 4, we briefly address the solitary wave 
solutions of the 6"^ GN model. Finally, Section 5 is devoted to conclusion. 


2. Derivation of the extended Green-Naghdi equations 
2.1. The extended GN system 

We hrst introduce the mean horizontal velocity component u = u{x,t) by 

1 n 

= (l)cc{x,y,t)dy, h = 1 + er/, (2.1) 

where h is the total depth of the fluid. The horizontal and vertical components of the 
surface velocity u and v are given respectively by 


yG)\y=eri^ ( 2 - 2 ) 

v{x,t) = (2-3) 

Multiplying (2.1) by h and differentiating the resultant expression by x and then using 
(1.1), (1.4), (2.2) and (2.3), we obtain the relation {hu)^ = erj^u — v/6^, or since erj^ = 

V = 5‘^{-{hu)^ + h^u). (2.4) 

Substitution of (2.4) into (1.2) yields the evolution equation for h: 

ht + e{hu)x = 0. (2.5) 

An advantage of choosing h and u as the dependent variables is that (2.5) becomes an 
exact equation without any approximation. 

Next, we differentiate (2.2) and (2.3) by x and t to obtain the relations 

(pxx T e(j)xy'f]x: Uf (f^xt T (2-6) 

W ^xy T ^4^yyVxj 4^yt T ^4^yyVtj (2-7) 

where the derivatives (pxx, (l>xy, (t>yy, (l>xt and (pyt are evaluated at y = erj. Similarly, 


{^(pt\y=erj^ X 4^xt T ^(pyiTjx. 


( 2 . 8 ) 
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Eliminating (jy^t and (j)yt with use of (2.6) and (2.7), (2.8) becomes 

(jj)t\y=er^X ^thx Vxht' (^■^) 

If we differentiate (1.3) by x, insert (2.4) and (2.9) in the resultant expression and use 
(2.5), we obtain the evolution equation for u\ 

Ut + Vthx + euUx + ehxUVx + rix = 0. (2.10) 

Using (2.5), Eq. (2.10) can be recast into the form 

[h{u + vhx)]t + e[h{u + vhx)u]x 

+e[hv{2hxUx + huxx) + h{u - u){ux + Vxhx)] + hr]x = 0 . ( 2 . 11 ) 

The system of equations (2.5) and (2.10) (or (2.11)) is equivalent to the basic Euler 
system (1.1)-(1.4). To obtain the extended GN equations, one needs to express the vari¬ 
ables u and V in (2.10) in terms of h and u. This is always possible as will be exemplihed 
below. Consequently, Eq. (2.10) can be written in the form ut = where 

Kn are polynomials of h and Unx,Unx,t, {unx = d'^u/dx'^,n = 0,1,2...). If one retains the 
terms up to order it yields the extended GN equation which is accurate to In 
accordance with this fact, we call the system of equations (2.5) and (2.10) (or (2.11)) with 
h and u being the dependent variables the extended GN system. 

2.2. The (5^ model 

Now, we derive the extended GN equation explicitly in the case of n = 2 by truncating 
the system of equations (2.5) and (2.11) at order 5^, which we call the model. To this 
end, we express the solution of the Laplace equation (1.1) subjected to the boundary 
condition (1.4) in the form of an inhnite series (see, for instance Whitham [15]) 

°° (y 1 'iSn a2n £ 

= 5 ——^ f2nx, f2nx = -K^, (2.12) 

where / = f{x,t) is the velocity potential at the fluid bottom y = —1. The expressions 
(2.1), (2.2) and (2.3) then become 

°° h‘2,n 

u = ( 2 . 13 ) 

n=0 ' ' 
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OO j ‘2iTi 

“ = E( 

n=0 ' ' 

(2.14) 

°o L2n-1 

n=l ^ ' 

(2.15) 


Retaining the terms of order S^, (2.13) gives 


U = fx- ^h^fxxx + ■^h'^fxxxxx + 0{5^). 
0 12U 


(2.16) 


The inverse relation in which is expressed in terms of h and u is achieved by the 
successive approximation starting from = u. This leads to 


^2 

fx u ~7rh Uxx T ^ 
6 


) --!^u 

l/t U^xxjxx -.^f^^XXXX 

oo 120 


+ 0 ( 6 ‘ 


(2.17) 


On the other hand, u from (2.14) is expanded as 


U = fx - ^h'^fxxx + ^h'^fxxxxx + O(S^). 

The last step is to introduce fx from (2.17) into (2.18), giving rise to the relation 


(2.18) 


5^ 


1 


^ ^ 3 ^ ^ ^ 9 ^ “1“ ^ ^xx'^xx )• (2.19) 

In view of (2.19), the expression of v from (2.4) can be written simply as 

(5R 


V = —5‘^hux - hfhxUxx + 0(5^). 

3 


( 2 . 20 ) 


The evolution equation for u follows by substituting (2.19) and (2.20) into (2.10) and 
using (2.5) to replace ht- After some manipulations, we arrive at the compact equation 
for u\ 


5^ 


45/i 


5 

ut + euux + ^ {h^{uxt + euuxx - 

[■{h {Uxxt T ^^^XXX he'Uj;'?/^;^;) 3 6/l A 0(^S 


( 2 . 21 ) 


If we multiply (2.21) by h and use (2.5), then we can put it in the conservation form 

/ /i2 \ ^2 

{hu)t + e \^hu^ + — J = — {h^{uxt + euuxx - 
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^4 

“^45 ^ + 0(5®). (2.22) 

The system of equations (2.5) and (2.21) (or (2.22)) is the extended version of the GN 
equation which is accurate to order 5^ and which retains the full nonlinearity. Recall 
that the amplitnde parameter e which characterizes the magnitnde of nonlinearity enters 
into h throngh the relation h = 1 + erj. What is meant by ”fnll nonlinearity” is that the 
parameter e is assnmed to be of order 1. 

If we use the approximation = 1 + enr]+ (e^) in (2.21) nnder the assnmption e 1, 
where n is a positive integer and ignore the e5^ and higher-order terms, it rednces to the 
classical Boussinesq system (Whitham [15]) 

hi di^hu^x b? T duUx ~\~ "djx ~^^xxt'! (2.23) 


and to the original GN system (Serre [1], Sn and Gardner [2] and Green and Naghdi [3]) 

^2 

ht + e{hu)x = ^, ut + euux + rix =— {h^{uxt + euuxx-eul)]^, (2.24) 

if we ignore the terms while preserving the fnll nonlinearity. 


Remark 1. Equation (2.21) involves the nonlinear terms of order e®. If one retains the 
and terms, then it rednces to a model eqnation presented in Madson and 
Schaffer [11] (see Eq. (3.11) therein). If one imposes a relation between the parameters e 
and 5, then one can derive the higher-order versions of the varions shallow-water model 
eqnations. For instance, the scaling e = 0(5^) 1 will lead to the higher-order KdV 

eqnation and to the higher-order GH eqnation nnder the scaling e = 0(5) 1. See 

Johnson [16] and Gonstantin and Lannes [17]. 


Remark 2. The linearization of the system of eqnations (2.5) and (2.21) abont the 
nniform state h = 1 and u = 0 yields the system of linear eqnations for r] and u 

52 _ 54 „ , , 

Vt T Ux 0, Ut -\- TJx ~Tr^xxt T ~7r'^xxxxt- (2.25) 

3 45 


Eliminating the variable rj from this system, one obtains the linear wave eqnation for u 

_ _ 52 _ 5 ^ _ 

Utt Uxx TT^xxtt T ~7^Uxxxxtt- (2.26) 

3 45 
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The dispersion relation for equation (2.26) is given by 


2 

3 45 

For the wave propagating to the right, the phase velocity Cp 
velocity Vg = doo/dk have the series expansions 


(2.27) 


oj/k and the group 


1 IQ 

c, = 1 - -{5kY + ^{5kf - i^{5kf + Omf), (2.28) 

1 iq 2c;q 

c. = 1 - ^{Skf + -{5kY - + 0{{Skf). (2.29) 

These expressions are compared with the exact dispersion relation u‘^ = k ianYikd/5, as 
well as the corresponding phase and group velocities which are given respectively by 


1 19 

Cp = 1 - -^{5kf + —{Skr - ^{Skf + 0 {{Skf), (2.30) 

1 19 

c, = 1 - -{Skf + ^{Skr - ^{Skf + 0{{Skf). (2.31) 

As expected, both velocities coincide up to order {kS^. Slight differences appear at order 

(M)F 


Remark 3. The system of equations (2.5) and (2.22) exhibits the following three con¬ 
servation laws: 


M= [h — l)dx^ 


' —OO 

/*oo 


H = 


hu^ + 

3 


hu dx, 
' 45 


h^u^x + ^ih-^y 


dx, 


(2.32) 

(2.33) 

(2.34) 


which represent the conservation of the mass, momentum and total energy, respectively. 
The first two conservation laws follow directly from (2.5) and (2.22) whereas the last one 
can be checked by means of a straightforward calculation using (2.5) and (2.21). A factor 
attached in front of the integral sign in (2.34) is only for convenience. 


3. Hamiltonian formnlation 



3.1. Hamiltonian structure 


Here, we show that the extended GN system (2.5) and (2.11) can be formulated as 
a Hamiltonian system by introducing an appropriate Lie-Poisson bracket. As pointed 
out by Zakharov [18], the basic Euler system of equations (1.1)-(1.4) conserves the total 
energy (i.e., kinetic plus potential energies) 


H = 


^2 roo r peri 


'-1 


K + ] dy + f] 


dx. 


(3.1) 


This fact can be conhrmed easily using (1.1)-(1.4). We substitute (2.12) into (3.1) and 
then perform the integration with respect to y to express H in terms of fx,fxx, ■■■ and 
h. Eliminating fx,fxx, ■■■ with the aid of the higher-order version of (2.17) yields a series 
expansion of H in (5^ as 


2 poo 


H = 


'y ^ fn{h, hx., hxxj ■■■, S Tj 


n=0 


dx = e^^5^^Hn, 


(3.2) 


n=0 


where fo = h and /n(> 1) are polynomials of h, hx, hxx, ■■■ (<5^ is simply a parameter). 
Specihcally, the hrst four Hn read 


H.Hr 


hu + - 1)' 


H. = 


dx, Hi = - h^uldx, H 2 = - — 


6 


90 


1890 


[^h^ulxx - 7h\hhx)xulx] dx. 


Let us now introduce the momentum density m = m{x, t) by 

5H 


em = 


5u' 


where the operator 5/5u is the variational derivative dehned by the relation 

d 


de 


H\u + etclL=o = 


.oo 

—— w dx, 
ou 


h u^^dx, 

(3.3) 

(3.4) 

(3.5) 


for arbitrary function w. In terms of m above, the Hamiltonian (or total energy) (3.2) 
can be put into the simple form 

roo 


H = 


eum + {h — 1)^ dx. 


(3.6) 
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Actually, it follows from (3.2) and (3.4) that 


m = 


n^nx Jnx' 


£^(-ir(/„ 

n=0 

Multiplying (3.7) by u and integrating by parts, we obtain 

OO ^CXD 

/ dx 6 ^ ^ hx-, hxx') 


(3.7) 


r*oo 


/ um ax = 
'-00 


(3.8) 


which, substituted in (3.2), gives (3.6). 

To proceed, we regard m and h as the dependent variables instead of u and h. Then, 
the variation of H in the former variables is found to be 


6H = 


eu Sm — 


OO p{n) 

n=0 s=0 


dfn _2 


dh 


}5h 


dx, (3.9) 


where p{n) is a nonnegative positive integer depending on n. In particular, p(0) = p(l) = 
0. To verify (3.9), we modify H from (3.6) as 


H = 


eum — -um H —(h — 1)^ 
2 2^ ^ 


dx, 


(3.10) 


and insert the expression of m from (3.7) into the middle term in the integrand, and then 
take the variations with respect to m, u and h. The coefficient of 6u vanishes by virtue of 
(3.7), resulting in (3.9). It readily follows from (3.9) that 

5H 


Sm 


= eu. 


SH 

Jh 


p{n) 


-VEE(-i)- 


<9/, 


dh 


n _2 


+ h-l. 


(3.11) 

(3.12) 


n=0 5=0 

Now, our main result is given by the following theorem. 


Theorem 1. The extended GN system (2.5) and (2.11) can he written in the following 
Hamiltonian form 


mt 

ht 


d^m + mdx hd^S fSH/Sm\ 

d^h 0 ) [SH/Sh J 


(3.13) 
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where the derivative = d/dx operates on all terms it multiplies to the right. Introduce 
the Lie-Poisson bracket between any pair of smooth functionals F and G 

r\5F^ ^ ^ JG SF^^dG 5F ^ dG] ^ ^ 

Then, the equations of motion form and h from (3.13) are expressed as 

mt = {m,H}, ht = {h,H}. (3.15) 


Proof. The proof of Theorem 1 relies on the following two relations: 

m = eh{u + vhx), (3.16) 

efivlflhxUx T huxx^ T h(u vfjiux T Vxhxff T hpx^ 

which we shall address hrst and then proceed to the proof of Theorem 1. 

We modify the Hamiltonian H from (3.1). First, use (1.1) to obtain 

dP 1 

~ (0x0)x + '^i4’y4’)y 

If we introduce this expression into (3.1) and perform the integration with respect to y 
under the boundary condition (1.4) as well as (2.2) and (2.3), we deduce 

/‘°° r/ V \ 2I 

H = — J yy-epxU + —jf) + 7] \^dx, 

where ip is the velocity potential at the free surface given by 

Ip = (p{x,er],t). (3.18) 
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Since this must hold for arbitrary m, we obtain an important relation 


m = ehipx, 


(3.19) 


which connects the variable ipx with the variable m. It follows by differentiating (3.18) 
by X and using (2.2) and (2.3) that 

'ipx = er], t) + e(f)y{x, erj, t)r]x = u + vh^, 

which, coupled with (3.19), yields (3.16). 

We substitute m from (3.16) into (3.17), divide the resultant expression by h and use 
(2.4). The relation thus obtained can be integrated by x to give 

1 6H 1 . . 1 


By taking the variational derivative of H from (3.1) with respect to h, we obtain 

.2 / 


(3.20) 


5H / 2 2 

«r=2 “ + 


V \ 5(1) 


dx. 


y=er} 


The following relation comes from the formula 5h{x, t)/5h{x', t) = 5(x—x'), where 5(x—x') 
is Dirac’s delta function: 


5(/)(x,y,t) 


5h(x', t) 


y=eri 


5)p{x, t) 
6h{x', t) 


v5{x — x'). 


We use this relation and (2.4) in the expression of 5H/5h to obtain 


5H 

Jh 


2?7 


-5)px 


. ^ H-+ e^{hu)xV + / hu^^dx. 

2 V 5^ e l_^ 5h 


Last, substituting the above expression into the left-hand side of (3.20), the relation to 
be proved reduces to 


ihu)xV -f f hu-ff-dx = hvUx — uu. 
. 5h 


(3.21) 


The integral on the left-hand side of (3.21) can be evaluated by using the relation ipx = 
m/eh from (3.19), which leads to 

r 

5h eh 
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showing that (3.21) holds identically. Hence, we establish (3.17). 

Now, Theorem 1 follows immediately from (3.16) and (3.17). First, we note that the 
introduction of (3.11) into (3.13) yields the system of evolution equations for m and h: 


rrit = —e{um)x — emu^ — h 



(3.22) 


ht = -e{hu)x. 


(3.23) 


Equation (3.23) is just Eq. (2.5). If we substitute (3.16) and (3.17) into (3.22), we see 
that the resultant equation is in agreement with Eq. (2.11). It is obvious that (3.15) with 
(3.14) is equivalent to (3.13). □ 


Remark 4. The bracket (3.14) has been used in Holm [19] to verify the Hamiltonian 
structure of two-dimensional shallow-water hydrodynamics with nonlinear dispersion, as 
well as that of the GN equation which has been detailed by Constantin [20]. It has 
a skew-symmetry {F,G} = —{G,F} and satishes the Jacobi identity {F, {G, H}} + 
{G, {H, F}} + {H, {F, G}} = 0 (Constantin [20]). 


Remark 5. If we retain the terms of order S^, then m and SH/Sh are found to be 


m = e 


6H e2 


^2 

hu — —{h^u. 
o 


x/x (h 


6h 2 

Equation (3.22) then becomes 


rrit = —euTTix — 2eUxm + —h 


^ - 4(h - 1) . 
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+ 5‘^h^ul - - ^(h - 1) 


(3.24) 

(3.25) 

(3.26) 


One can check by a direct computation that (3.26) indeed reproduces (2.22). 


Remark 6. In addition to the conservation of mass, momentum and total energy which 
are given respectively by (2.32), (2.33) and (2.34), the extended GN system exhibits the 
fourth conservation law L = J^^m/hdx. This immediately follows from the equation 
{m/h)t = —{eum/h + SH/6h)x which is derived from the system of equations (3.22) and 
(3.23). In view of (3.16), the density of L is written as m/h = e{u + vh^)- Since the 
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unit tangent vector t to the free surface is given by t = (1/ a /1 + hx/ a /1 + hi) and 
the surface velocity by u = {u,v), the density m/h is equal to the tangential component 
= u ■ t of the surface velocity multiplied by a factor eA/TT~^. It turns out that L 
is represented by the line integral of the tangent velocity along the free surface. To be 
more specihc, L = e J^Utds, where ds = \/l + h^dx is the line element of the curve C 
representing the free surface. Using (3.19), this integral can be evaluated explicitly as 
L = e('^_i_ — with = hma;_j.±oot/’(a^U)) showing that the conserved quantity L is 
determined only by the limiting values 'ip± of the velocity potential at the free surface. 
Note also that the quantity mdx is conserved. However, it is not independent of the 
momentum P given by (2.33). Actually, integrating the relation (3.7) with /q = h, the 
above conserved quantity reduces to e hu dx = eP. 

3.2. Relation to Zakharov’s Hamiltonian formulation 


The water wave problem posed by (1.1)-(1.4) permits a Hamiltonian formulation. 
Following Zakharov [18], the equations of motion for surface gravity waves can be written 
in terms of the variables h and tjjx as 


ht = -- 


'ipxt = 


6H 


5jr\ 
Jh ). 


Note that the variable ipx is used instead of fj. This choice 
analysis. If we define the bracket between any pair of smooth 


(3.27) 

(3.28) 

is suitable for the present 
functionals F and G by 


{F,G} 



_Sh XdfjxJ 


then, Eqs. (3.27) and (3.28) are expressed as 


( dF\ 5a 
\5fJx) X 


(3.29) 


ht = {h, 77}, fjxt = {fJx, H}- 


(3.30) 


Under the above setting, the following theorem is established: 


Theorem 2. The system of equations (3.13) is equivalent to Zakharov’s system of equa¬ 
tions (3.27) and (3.28). 
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Proof. First, we change the variables h and to h and m and rewrite Eqs. (3.27) 
and (3.28) in terms of the latter variables. To this end, we evaluate the variation of the 
Hamiltonian H in two alternative ways: 


poo 

Jh 

SH 


poo 

'5H 


- 

/ 



5ipx 

dx = 


5h+ — 

5m 

/ —oo 

5h 

V)^:fixed "17 a: 

/iifixed 

J —OO 

5h 

m:fixed 

/iifixed 


(3.31) 

We compute the variation of m from (3.19) as 5m = et/jxSh + ehdipx and subsitute this 
expression into the right-hand side of (3.31) and then compare the coefficient of 5h and 
5'il)x on both sides. It follows that 


5H 

Ih 




6H 

Jh 


m: fixed 


6H 

5m 


'4^X1 

h: fixed 


(3.32) 


5H 


5'ipx 


fixed 




fixed 


Taking into account (3.11) and (3.19), the above relations become 


5H 

Jh 


fixed 


5H 

Jh 


mifixed 


eum 


(3.33) 


(3.34) 


5H 


54^x 


= e^hu. 


(3.35) 


/i; fixed 


These formulas enable us to transform the Zakharov system to the extended GN 
system. Actually, Eq. (3.27) with (3.35) yields Eq. (3.23). On the other hand, introducing 
(3.19) and (3.34) into (3.28), we obtain 


m\ mt mht f 5H\ /um\ 
h)t h X ^ h )x 


(3.36) 


where we have used the notation (5i7/(5h)|m:fixed = 5H/5h for simplicity. If we replace ht 
in the middle expression of (3.36) by the right-hand side of Eq. (3.23), then the equation 
reduces to Eq. (3.22) after multiplying h on both sides. Last, we use (3.32) and (3.33) 
with H being replaced by F and G, respectively and (3.19) for ipx to rewrite the bracket 
(3.29) in terms of h and m. The resultant expression becomes 


{F,G} 



5J^ m5F\ 
5h^ h 5m) 



(h—'] (— 

\ 5m) x\^h h 5m) 


dx. 


(3.37) 
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Integrating the second term of (3.37) by parts and rearranging the integrand, we conhrm 
that (3.37) transforms to the bracket dehned by (3.14). □ 

4. Solitary wave solutions 

We seek the solitary wave solutions of the extended GN system (2.5) and (2.22) in 
the form of the traveling wave h = h{^),u = u{^), = x — ct). We impose the boundary 

conditions h — )■ 1, h' —)■ 0 and u,u',u",u"',u'"' — )■ 0 as |a;| — )■ oo, where c(> 0) is the 
velocity of the solitary wave and primes refer to differentiation with respect to If we 
introduce these forms into (2.5) and (2.22) and integrate each equation once with respect 
to ^ while taking into account the boundary conditions, we deduce 


-chu + € hu^ + 


r 

+ 45 


—c{h — 1) + ehu = 0, 
h^-l\ 6^ 


2e2 


= —h^{—cu" + euu" — e{u'y) 
3 


[h^{—cu"' + euu'" — 5eu'u")y — 3eh^{u"y 
We use (4.1) to express u in terms of h as 



(4.1) 


(4.2) 


(4.3) 


Substituting (4.3) into (4.2), multiplying the resultant equation by h'/h? and then inte¬ 
grating, we obtain the third-order nonlinear differential equation for h\ 


— (c^ -I- 2)h^ -I- (2c^ -I- l)h — 

= ^ Wh'h'" - h\h"f + 2h{h'fh" - I2{h'f] . (4.4) 

3 45 

If we neglect the <5^ terms on the right-hand side of (4.4), then it reduces to the 
corresponding equation for the GN equation. Unlike the Boussinesq system (2.23), it 
exhibits the analytical solitary wave solution of the following form which has been obtained 
for the hrst time by Serre [1] and later by Su and Gardner [2]: 


h 


1 -1- (c^ — 1) sech' 


:V/3F^ 


2c5 


e, (ol). 


(4.5) 
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An inspection reveals, however that Eq. (4.4) would not have analytical solutions and 
hence numerical analysis may be necessary to extract the characteristics associated with 
the solitary waves such as the surface prohle and the velocity-amplitude relation. Instead, 
we perform a singular perturbation analysis. 

Note hrst that the parameter S in Eq. (4.4) disappears by means of the scaling ^ —)■ S^. 
Furthermore, the form of the solitary wave solution (4.5) suggests that an appropriate 
scaling of the variable ^ would be b(, where C = ^/S and the unknown parameter b plays 
the role of the wavenumber which is determined in the course of the pertubation analysis. 
In accordance with this observation, we expand h and in powers of e as 


h{b() = 1 + ehi + e^/i2 + O(e^), 

(4.6) 

= 1 -|- eci -f- e^C2 -|- e^C3 -|- O(e^), 

(4.7) 

= ebi + e%2 + e% + 0{e^). 

(4.8) 


We employ the conditions hi(0) = 1 and hj{0) = 0,j > 2 which dehne the amplitude e 
by /i(0) = 1 -|- e. The method for constructing solitary wave solutions is well-known. See, 
for instance Fenton [21]. Hence, we summarize the result. 

The solution to the boundary value problem includes the secular term like ( sech^b^ tanh b^. 
To ensure the uniform validity of the expansion (4.6), one must eliminate the secular terms 
by choosing the coefficients bj and Cj (j = 1,2,...) appropriately. The solitary wave solu¬ 
tion thus constructed can be written as 

Q 5 / 14Q \ 

h{b() = l+€sech.‘^b(—-€‘^sech.‘^b(ta.nh^b(+-€^ f 1 — —sech.^b(j sedi^bCtanh^ b(+0{€'^), 

(4.9) 

= 1 + e - - ie’ + 0(/), (4.10) 

The above solution coincides with the third-order solution obtained in Fenton [21] and 
Grimshaw [22] up to order e^. The discrepancy appears at order e^. The reason is obvious 
since the latter takes into account the <5® terms which have been neglected in our model 
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equation. To be more specific, the scaling e = 0(5^) is assumed in Fenton [21] and 
Grimshaw [22] so that the 5® terms turn out to be comparable with the terms. 

Remark 7. Introduce the quantity Q by 



which has been shown to be conserved by the extended GN system (See Remark 6). For 
the model, the solitary wave solntion h{^) and m(.^) is a critical point of the fnnctional 
H — cQ with respect to m and h. Specifically, 


5m ^ 5m ’ 5h ^ 5h 


(4.13) 


This fact can be confirmed by a direct compntation nsing (3.11), (3.24), (3.25) and (4.4). 
Thus, a variational characterization is possible for the solitary wave solntion of the 5‘^ 
model, even if it does not exhaibit exact solitary wave solutions. It will be an important 
issne to establish this assertion for the general 5^” (n > 3) model. In the case of the GN 
eqnation, the similar result has been pointed out and used to prove the linear stability of 
the solitary wave solntion (4.5) (Li [23]). 


5. Conclusion 

We have presented a novel method for extending the GN eqnation to the general sys¬ 
tem which incorporates the arbitrary higher-order dispersive terms while preserving the 
fnll nonlinearity. As an illustrative example, we have derived a model equation which is 
accnrate to order for which the perturbation analysis reprodnces the known solitary 
wave solntions at the level of order e^. Nevertheless, the nnmerical analysis will be neces¬ 
sary to extract the fnlly nonlinear nature of the model equation, as has been performed 
for the GN eqnation by Li et al [24] and Mitsotakis et al [25]. We have shown that our 
extended system permits the same Hamiltonian strnctnre as that of the GN eqnation. 
In the process of establishing this fact, the energy integral of the basic Euler system has 
been used efficiently. We have also verified that Zakharov’s Hamiltonian formulation of 
the cnrrent water wave problem is eqnivalent to that of the extended GN eqnations. This 
has been accomplished by rewriting Zakharov’s eqnations in terms of the total depth of 
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the fluid and momentum density. The extension of our model equations to the more gen¬ 
eral setting with an uneven bottom topography and its three-dimensional generalization 
will be done along with the same procedure as that developed here. In any case, there 
are many things to be resolved by future study. 
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